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Preface 


A  preface  to  the  entire  Handbook  of  Supersonic  Aerodynamics 
appears  in  Volume  1  and  includes  a  brief  history  of  the  project.  As 
stated  in  Volume  1,  "The  primary  criterion  used  in  selecting  material 
for  the  Handbook  is  its  expected  usefulness  to  designers  of  supersonic 
vehicles.  Thus  a  collection  of  data  directly  useful  in  the  design  of 
supersonic  vehicles,  results  of  the  more  significant  experiments,  and 
outlines  of  basic  theory  are  included  ...." 

The  present  edition  of  the  Handbook,  printed  and  distributed 
by  the  Bureau  of  Ordnance,  is  being  published  in  separate  sections  as 
material  becomes  available.  The  contents  of  the  entire  work  is  given 
on  the  back  of  the  title  page.  It  may  be  noted,  the  majority  of  sec¬ 
tions  of  the  Handbook  are  being  prepared  by  individual  authors  for  the 
Aoplied  Physics  Laboratory.  Those  sections  are  noted  for  which  an 
early  publication  date  is  expected.  The  selection  of  material,  edit¬ 
ing,  and  technical  review  of  the  Handbook  of  Supersonic  Aerodynamics 
continue  to  be  carried  out  by  an  editor  and  a  technical  reviewing  com¬ 
mittee  at  the  Applied  Physics  Laboratory. 

Volume  3  of  the  Handbook  of  Supersonic  Aerodynamics  contains 
three  closely  related  sections,  "Two-Dimensional  Airfoils"  (Section  6), 
"Three-Dimensional  Airfoils"  (Section  7),  and  "Bodies  of  Revolution" 
(Section  8).  A  fourth  related  section  on  wing-body  interference 
(Section  9)  is  placed  in  Volume  4  with  sections  on  stability  and  con¬ 
trol  because  of  its  importance  to  this  later  subject. 


i- 

□ 

□ 


The  staff  of  the  Aerodynamics  Handbook  Project  at  the  Applied 
Physics  Laboratory  gives  grateful  acknowledgment  to  the  Aeromechanics 
Division  of  the  Defense  Research  Laboratory  of  The  University  of  Texas 
for  the  large  amount  of  computational  work  carried  out  by  them  for  the 
Sections  6  and  7  of  the  Handbook. 

The  numbering  system  of  Volume  3  is  essentially  the  same  as 
that  used  in  the  last  preceding  volumes. 


August  1957 


NAVORD  REPORT  1488  (Volume  3) 
ii 


V 

V  \ 


Prefac 


A 


T 


Agencies  and  individuals  interested  in  the  aeronautical  sci¬ 
ences  are  invited  to  submit  and  to  recommend  material  for  inclusion 
in  the  Handbook;  full  credit  will  be  given  for  all  such  material  used. 
Regarding  the  selection  of  material  and  the  preparation  of  the  volumes 
in  the  Handbook  Series,  the  Applied  Physics  Laboratory  earnestly  soli¬ 
cits  constructive  criticisms  and  suggestions.  Correspondence  relating 
to  the  editing  of  the  Handbook  should  be  directed  to 


Editor  and  Supervisor,  Aerodynamics  Handbook  Project 

Applied  Physics  Laboratory 

The  Johns  Hopkins  University 

8621  Georgia  Avenue 

Silver  Spring,  Maryland 


Communications  concerning  distribution  of  the  Handbook  should 
be  directed  to 


Chief,  Bureau  of  Ordnance 
Department  of  the  Navy 
Washington  25,  D.  C. 
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SECTION  7  -  THREE-DIMENSIONAL  AIRFOILS 
Primary  Symbols 

velocity  of  sound 
b2 

aspect  ratio:  g— 


wing  span  (twice  the 
distance  from  the 
tip  to  the  root 
chord) : 


wing  chord 


T>r 


5<r 


— ^ — 

:,-T 


drag  coefficient: 


Drag  Force 


lift  coefficient: 


Lift  Force 


lift-curve  slope: 


pitching  moment  coefficient: 

pitching  moment  about  the 
„  _  apex  of  the  leading  edge 


3PV  Scr 

(where  a  positive  moment  is  defined 
as  a  moment  tending  to  tilt  leading 
edge  upward) 

dC 

moment-curve  slope: 


P  "  P. 


pressure  coefficient: 


O  |c* 
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E 


(  &  cot  A  ) 


g 

M 

P 

S 

t 


complete  elliptical  integral  of  the  second 
kind,  of  modulus  Vi  -  (  0  cot  A  ) 

source  intensity  per  unit  area 

V 

Mach  number:  — 
a 

pressure 

wing  planform  area  (two  panels) 
thickness  of  wing  profile 


thickness  ratio 


u,  v,  w 


components  of  perturbation  velocity  in 
x-,  y-,  z-  directions 


V 


velocity  of  the  undisturbed  stream 


x,  y,  z 


rectangular  coordinates 


a 

P 

6 


A 


A 


T 


X 


moment  arm  of  wing  elemental  lifting  area 
angle  of  attack 

|/m2  -  1 

tip  rake  angle  (angle  which  tip  makes  with  the 
root  chord  direction;  defined  as  positive 
when  the  tip  is  raked  inboard) 


leading  edge  sweepback  (angle  which  the  lead¬ 
ing  edge  makes  with  normal  to  root  chord 
direction) 

trailing  edge  sweepback  (angle  which  the  trail¬ 
ing  edge  makes  with  normal  to  root  chord 
direction) 

taper  ratio  (tip  chord  divided  by  root  chord 
for  unraked  wings  with  straight  leading 
and  trailing  edges) 


Mach  angle: 


sin 


-1 


I 

M 
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Introduction 


700.1 


SECTION  7  -  THREE-DIMENSIONAL  AIRFOILS 


This  section  of  the  Handbook  of  Supersonic  Aerodynamics  was 
prepared  at  the  Applied  Physics  Laboratory  of  The  Johns  Hopkins  Uni¬ 
versity.  Many  of  the  formulae  and  graphs  for  lift  and  moment  charac¬ 
teristics  of  three-dimensional  airfoils  were  especially  prepared  for 
this  section  by  the  Defense  Research  Laboratory  of  The  University  of 
Texas. 


700  Introduction 


This  section  of  the  Handbook  presents  the  aerodynamic 
characteristics  of  a  limited  class  of  finite  wings  in  a  convenient 
form  for  use  in  design  calculations.  The  treatment  of  finite  wings 
is  severely  limited  to  uncambered  wings  with-  zero  angle  of  attack 
for  which  thickness  ratio  and  position  of  maximum  thickness  are  held 
at  a  constant  value.  Section  6  of  the  Handbook  (Two-Dimensional  Air¬ 
foils)  may  be  used  for  guidance  in  estimating  viscosity  effects  and 
the  effect  of  changing  the  angle  of  attack,  thickness  ratio,  camber, 
and  position  of  maximum  thickness. 

Section  7  treats  wings  alone  and  ignores  the  fact  that,  in 
practical  application,  wings  are  attached  to  bodies.  Bodies  alone 
are  discussed  in  Section  8  of  the  Handbook  (Bodies  of  Revolution), 
and  interference  effects  of  the  wing-body  combination  are  treated  in 
detail  in  Section  9  of  the  Handbook  (Mutual  Interference  Phenomena) . 

Section  7  also  includes  a  summary  of  the  underlying  theory 
with  particular  emphasis  on  the  limiting  conditions  which  affect  the 
range  of  applicability  of  the  data  presented. 


700.1  Scope  of  Contents 

Theoretical  airfoil  characteristics  for  lift,  pitching 
moment,  and  pressure  drag  are  computed  by  the  method  of  supersonic 
source  distributions.  They  are  presented  in  graphical  form  using 
non-dimensional  parameters  which  define  the  aerodynamic  and  geometric 
properties  of  wings  in  supersonic  flow.  The  lift  and  moment  analysis 
is  discussed  first  because  much  of  it  is  also  a  basic  consideration 
in  the  drag  analysis.  The  calculations  at  zero  angle  of  attack  in¬ 
volve  linearized  theory  only.  The  reader  is  referred  to  Section  6 
of  the  Handbook  (Two-Dimensional  Airfoils)  for  a  comparison  of  lin¬ 
earized  theory  and  higher  order  theories. 

Throughout  the  discussion  theoretical  analyses  are  reduced 
to  a  minimum  consistent  with  intelligent  use  of  the  data,  and  a 
reference  is  made  to  the  original  sources  containing  the  more  de¬ 
tailed  analyses.  Sample  calculations  are  included  to  demonstrate 
the  practical  use  of  the  graphs  in  determining  the  aerodynamic  char¬ 
acteristics  of  several  wing  configurations  at  specified  flight  con¬ 
ditions  . 
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A  brief  comparison  of  theory  and  experiment  is  included  in 
order  to  provide  a  convenient  basis  for  estimating  the  accuracy  and 
reliability  of  theoretical  predictions  from  the  linearized  theory. 
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Resume  of  Basic  Theory 


701.2 


701  Resume  of  Basic  Theory 

A  resume  of  basic  theory  is  sketched  out  here  to  assist  the 
reader  in  his  understanding  and  interpretation  of  the  material  pre¬ 
sented  in  graphical  form.  This  discussion  is  limited  to  a  considera¬ 
tion  of  a  supersonic  source  solution.  The  more  complex  questions, 

e.g.,  boundary  conditions,  conical  flow  solutions,  and  so  forth,  are 
referenced  where  appropriate. 


701.1  Basic  Flow  Assumpt ions 

The  subsequent  analysis  is  based  on  the  assumption  of  a  fluid 
prescribed  by  the  following  characteristics: 

a.  Compressible  fluid 

b.  Inviscid  fluid,  without  heat  transfer 

c.  Absence  of  all  external  force  fields 

d.  Steady  motion,  i.e.,  flow  quantities  are  independent  of 

time 

e.  Irrotational  flow,  i.e.,  the  flow  field  is  free  of 

vorticity  (see  item  g  below) 

f.  Shock  waves  are  infinitely  weak 

g.  Adiabatic  processes,  i.e.,  no  heat  transfer  at  the 

boundaries 


701.2  General  Equation  of  Motion 

The  equation  of  motion  governing  the  flow  under  these  condi¬ 
tions  is  expressed  in  the  total  velocity  potential  form  (Ref.  1, 
p.  2.17,  Eq.  2.31)  as  follows: 


(701.2-1) 


where  is  the  total  velocity  potential  defined  by 


and 


etc 
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The  primed  symbols  represent  the  total  velocity  components  as  con¬ 
trasted  with  perturbation  velocities  which  are  designated  by  symbols 
without  primes. 


Here  a 


the  velocity  of  propagation  of  an  infini¬ 


tesimal  pressure  wave  at  constant  entropy  (sonic  velocity). 


701.21  Linearized  Equation  of  Motion 

Solutions  for  this  second  order  non-linear  differential 
equation  (Eq.  701.2-1)  have  been  found  for  only  a  limited  number  of 
cases  and,  therefore,  some  simplification  is  necessary  if  any  practi¬ 
cal  applications  are  to  be  made.  In  applications  to  the  flow  past 
finite  wings,  a  large  range  of  important  wing  configurations  meet  the 
requirements  of  the  linearization  of  this  equation.  This  lineariza¬ 
tion  is  a  consequence  of  using  the  small  perturbation  theory  in  which 
it  is  assumed  that  the  disturbances  produced  by  a  body  in  a  fluid  are 
very  small  compared  to  the  velocity  of  the  undisturbed  flow.  The 
wings  of  aircraft  flying  at  supersonic  speeds  are  usually  sufficiently 
thin  so  that,  at  small  angles  of  attack,  the  velocity  disturbances  can 
be  assumed  to  satisfy  the  requirements  of  the  small  perturbation  theory 
at  all  points  in  the  flow,  except  in  the  immediate  vicinity  of  the  sub¬ 
sonic  leading  edge. 

The  linearization  of  Eq.  701.2-1  by  the  theory  of  small  per¬ 
turbations  is  accomplished  as  follows: 

a.  The  total  velocity  components  u' ,  v’,  wr  in  the  undis¬ 

turbed  stream  are: 

u*  =  V  (constant) 

v'  =  0 

w'  =  0 

b.  The  corresponding  total  velocity  components  in  the 

presence  of  the  disturbance  are: 

u'  =  V  +  u 

v'  =  V 

w'  =  w 

where  u,  v,  and  w,  the  disturbance  velocities,  are 

limited  by  the  small  perturbation  theory  to 

u  «  V,  v  «  V,  w «  V. 

Since  the  flow  external  to  the  disturbing  body  is  assumed  irrotational, 
a  perturbation  velocity  potential,  # ,  is  introduced  defined  by 

=  t(x,  v  =  tfy,  w  =  flfz. 


u 
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Substituting  in  the  general  potential  Eq.  701.2-1,  and  neglecting 
terms  of  second  order  and  higher  in  the  perturbation  velocities  and 
their  first  derivatives,  we  arrive  at  the  following  linearized  form 
of  the  equation  of  motion: 


Moo2  ~  1 


0....  -  "  0  =  0 


xx 


yy 


zz 


(701.21-1) 


v 

where  M  <*,  =  - —  is  the  Mach  number  of  the  undisturbed  flow.  This 

a  OO. 

linearization  leads  to  solutions  of  the  same  degree  of  accuracy  as 
that  given  by  the  first  order  theory  for  two-dimensional  airfoils 
which  was  treated  in  Section  6  of  the  Handbook. 


701.3  Solut ion  of  Linearized  Equation  by  the  Method  of  Su  jnic 
Source  Distributions 

An  elementary  solution  of  the  linearized  form  of  the  equa¬ 
tion  of  motion  (Eq.  701.21-1)  is  defined  by  the  following  function 
(see  Puckett,  Ref.  2,  Eq.  3): 


<A  (x,  y,  z)  = 


-C 


V(x  -C)2  "  P2  [(y  -  v)2  +  (z  -C)2] 

(701.3-1) 


where 


/32  =  mJ  -  1  . 


The  function  defined  by  Eq.  701.3-1  represents  the  potential 
at  the  point  (x,  y,  z)  due  to  a  point  source  of  strength  proportional 
to  C  and  located  at  the  point  P  (£  ,  77,  £  )  (see  Heaslet  and  Lomax, 

Ref.  3,  p.  149).  This  function  becomes  imaginary  outside  the  cone 
shown  in  the  figure  on  the  next  page  and  described  by  the  equation, 

(x  -  £)2  -  P2  [(y  '  t?)2  +  (z  -  C)2]  =  0.  (701.3-2) 

This  cone  has  its  vertex  at  the  point  P  (  £  ,  T),  t,  )  and  has  a 
semi-vertex  angle  of 


\i  =  tan-1  ^  =  sin  1  i  . 


It  corresponds  to  the  Mach  cone  formed  by  the  stationary  wave  front 
produced  in  a  flow  at  Mach  number  M  <*,  by  a  disturbance  at  P  (  £  ,  r\,  £  ) . 


4890 JO  0-59-2 
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A  disturbance  at  P  (  (  ,  T),  5  )  cannot  be  felt  outside  the  cone  down¬ 
stream  from  this  point,  but  is  felt  within  the  cone,  thereby  defining 
zones  of  inaction  and  action  respectively.  The  zone  of  action  is 
noted  as  the  Zone  of  Influence  in  the  figure  above.  Likewise,  condi¬ 
tions  at  P  (  (  ,  Vt  t  )  depend  only  on  events  occurring  within  the  cone 
upstream  of  the  point.  This  zone  is  known  as  the  Zone  of  Dependence. 

Solutions  of  more  general  applicability  are  obtained  by  as¬ 
suming  distributions  of  sources  of  varying  strength.  This  is  permis¬ 
sible  since  the  differential  equation  is  linear  and  the  principle  of 
superposition  applies.  Puckett  (Ref.  2)  has  applied  the  distributed 
source  solution  to  the  evaluation  of  the  flow  over  a  thin  triangular 
planform  wing.  The  analysis  is  summarized  here  as  a  means  of  illus¬ 
trating  the  application  of  the  small  perturbation  theory  to  the 
determination  of  wing  characteristics  at  supersonic  speeds.  The  thin 
wings  under  consideration  are  assumed  to  lie  essentially  in  the  xy-plane 
and  the  sources  are  distributed  in  that  plane.  Letting  g  represent  the 
source  strength  per  unit  area,  the  differential  element  of  the  pertur¬ 
bation  velocity  potential,  d0,  is  obtained  by  writing  (g  d£  dT?)  in 
place  of  C  in  Eq.  701.3-1.  Therefore, 

d£f  =  —  ~S  -  ,  (701.3-3) 

(x  "  £  )2  "  Z32  [  (y  -  Tj)2  +  z2] 


where 


g  =  g(  i  ,  *?) 
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701.31  Determination  of  Velocity  Potential  for  Thin  Wings 

The  velocity  potential  at  any  point  on  the  wing  surface  is 
obtained  by  the  integration  of  Eq.  701.3-3  over  the  wing  area  in¬ 
closed  within  the  Zone  of  Influence  for  that  point.  Puckett  (Ref.  2) 
demonstrated  that  the  source  intensity,  g,  at  any  point  is  dependent 
only  upon  local  conditions,  and  is  given  by 


g  =  J  .  (701.31-1) 

The  boundary  condition  simply  requires  that  the  local  velocity  be 
parallel  to  the  wing  surface  at  any  point,  that  is, 

u >x  =  v  ’  (701.31-2) 


where  u>x  is  the  slope  of  the  wing  surface  relative  to  the  x-direction 

at  a  given  point.  Since  these  wings  are  assumed  to  lie  in  the  xy-plane, 
it  is  sufficient  to  satisfy  the  boundary  condition  in  that  plane 
(z  =  0).  The  velocity  potential  at  the  point  (x,  y)  becomes,  after 
substitution  from  Eqs.  701.31-1  and  701.31-2  in  Eq.  701.3-3, 


tf(x,  y) 


d  £  d7} 


-/8  2  (y  -  r/) 2 


(701.31-3) 


where  AS  is  the  area  of  source  distribution  intercepted  by  the  fore¬ 
cone  from  the  point  (x,  y)  and  thus  lies  within  the  Zone  of  Influence. 


701.32  Perturbation  Velocity  and  Pressure  Coefficient 

The  horizontal  perturbation  velocity  in  the  x-direction,  $x, 
is  obtained  from  Eq.  701.31-3.  The  pressure  coefficient,  defined  as 


c 


P 


2 


is 


c 

P 


-2u 
V  ’ 


(701.32-1) 


(701.32-2) 


as  obtained  from  the  linearized  theory 
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701.4  Infinite  Triangular  Wing 

The  infinite  triangular  wing  is  used  as  a  basic  configuration 
in  supersonic  flow  analysis  somewhat  analogously  to  the  infinite  aspect 
ratio  wing  in  subsonic  flow  analysis.  The  properties  of  finite  wings 
in  supersonic  flow  are  frequently  obtained  by  the  superposition  of  tri¬ 
angular  wing  components. 


701.41  Supersonic  and  Subsonic  Leading  Edges 

In  supersonic  wing  theory  the  treatment  of  supersonic  flow 
over  wing  surfaces  is  divided  into  two  main  divisions.  These  divisions 
are  referred  to  as  the  supersonic  and  subsonic  leading  edge  conditions 
in  which  the  velocity  component  normal  to  the  leading  edge  is  super¬ 
sonic  and  subsonic  respectively.  These  conditions  are  also  identified 
by  the  position  of  the  leading  edge  with  respect  to  the  Mach  line  from 
a  point  on  the  leading  edge  (as  shown  in  the  figure  below). 


lo)  Supersonic  Leading  Edge  lb)  Subsonic  Leading  Edge 


With  supersonic  leading  edges,  pressure  changes  due  to  dis¬ 
turbances  on  one  surface  cannot  be  transmitted  forward  and  thence 
around  the  leading  edge  to  modify  the  flow  over  the  other  surface. 

The  flow  over  upper  and  lower  surfaces,  therefore,  can  be  treated 
independently.  On  the  other  hand,  a  subsonic  leading  edge  condition 
permits  the  transmission  of  pressure  disturbances  from  one  surface  to 
the  other  and  these  interferences  must  be  taken  into  consideration  in 
the  analysis  of  the  flow  over  the  wing  surfaces. 
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701.42  Horizontal  Perturbation  Velocity  for  a  Triangular  Wing 

Puckett's  analysis  (Ref.  2)  treats  the  infinite  triangular 
wing  with  a  supersonic  leading  edge.  The  Mach  line  from  the  wing  apex 
in  this  instance  divides  the  wing  surface  into  two  regions,  and  the 
integration  (indicated  by  Eq.  701.31-3)  is  accomplished  separately  for 
each  of  these  regions.  The  horizontal  perturbation  velocity,  u,  is 
constant  in  the  region  ahead  of  the  Mach  line  (Region  Aj,  Fig.  (a)  on 

the  preceding  page),  and  is  given  by  the  following  equation  (Puckett, 
Ref.  2,  Eq.  34): 


-u)x  ( P  cot  A  ) 

P  /( P  cot  A  ) 2  -  1 


(701.42-1) 


In  the  region  behind  the  Mach  line  (Region  Ag,  Fig.  (a)) 
tion  velocity  is  given  by  the  equation  (Ref.  2,  Eq.  35): 


the  perturba- 


“2wx  (0  cot  A  )  /  i  -£2(*)2 

-  - ■— . - . .  cos  / - 

npV(P  cot  A  )2  -  1  V  (P  cot  A  )2  -  p2($)2 


(701.42-2) 

Stewart,  using  the  method  of  conformal  transformation,  deter¬ 
mined  the  flow  characteristics  for  an  infinite  triangular  wing  with 
subsonic  leading  edge  and  his  results  (Ref.  5,  Eq.  45),  rewritten  in 
the  present  notation,  give 


u 

V 


-u>x  (P  cot  A  )2 
P  P  cot  A  )2  -/32(*)2 


1 


7 


( P  cot  A  ) ' 


(701.42-3) 


where  *  -  ( P  cot  A  )2/  is  the  complete  elliptic  integral  of  the 

second  kind,  of  modulus  -v/l  -  ( P  cot  A  ) 2.  Equations  701.42-2  and 
701.42-3  show  that  the  perturbation  velocity,  u,  and  hence  the  pressure 

coefficient,  c^  =  — 2u,  remain  constant  along  radial  lines  drawn  from 

the  apex,  which,  according  to  Busemann  (Ref.  6),  is  characteristic  of 
conical  flow  phenomena. 
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701.5  Analysis  of  Finite  Wing  Planf orms 

The  aerodynamic  characteristics  of  thin,  flat,  sweptback 
finite  wings  are  di .ormined  by  suitable  modification  and  application 
of  the  results  for  infinite  flat  triangular  wings.  When  the  boundary 
of  a  finite  configuration  crosses  the  basic  triangular  wing,  a  correc¬ 
tion  must  be  determined  which  in  effect  cancels  the  pressure  field  of 
the  infinite  triangular  wing  beyond  the  boundary  of  the  finite  configu¬ 
ration  without  altering  the  specified  boundary  conditions.  In  some 
instances,  e.g.,  for  a  triangular  planform,  it  is  sufficient  to  locate 
the  trailing  edge  boundary  and  disregard  the  downstream  portion  of  the 
pressure  field  of  the  wing,  since  pressure  disturbances  cannot  be  trans 
mitted  upstream.  An  example  of  this  type  of  analysis  is  given  in  the 
computation  of  the  drag  of  a  triangular  planform  with  a  double  wedge 
cross  section  (see  Subsection  702.22). 

Other  planforms  having  supersonic  trailing  edges  require  more 
elaborate  and  detailed  pressure-field  cancelling  techniques  to  account 
for  the  effect  of  disturbances  radiating  from  planform  discontinuities 
such  as  the  leading  edge  of  a  tip  chord.  The  flow  about  wings  moving 
at  supersonic  speeds  with  subsonic  trailing  edges  must  satisfy  the  well 
known  Kutta  condition  which  requires  the  velocity  to  be  finite  at  the 
trailing  edge. 

The  analysis  and  data  presented  herein  are  limited  to  wings 
having  sweptback  leading  edges  only.  Data  for  wings  with  sweptforward 
leading  edges  can  be  determined  by  means  of  the  reversibility  theorem, 
discussed  later  (see  Subsection  702.3),  from  data  for  equivalent  wings 
having  sweptback  leading  edges. 

The  source  distribution  methods  cited  in  Puckett  (Ref.  2) 
and  Evvard  (Ref.  7)  offer  the  most  convenient  means  of  calculating  the 
aerodynamic  characteristics  for  planforms  with  supersonic  leading  and 
trailing  edges.  The  calculation  of  the  flow  properties  for  wings  with 
subsonic  leading  edges  is  accomplished  through  the  use  of  solutions 
presented  in  Stewart  (Ref.  5)  and  Lagerstrom  (Ref.  8)  for  the  infinite 
triangular  wing  plus  suitable  pressure  cancelling  techniques.  Two 
such  methods  or  techniques  appear  to  be  generally  used.  One  of  these 
methods,  discussed  by  Mirels  (Ref.  9),  uses  supersonic  doublets  and 
is  applicable  to  curved  boundaries.  The  second  method,  discussed  by 
Cohen  (Ref.  10),  uses  the  superposition  of  conical  flow  fields  origi¬ 
nally  introduced  by  Busemann  (Ref.  6)  .  The  details  of  these  pressure- 
field  cancellation  techniques  are,  however,  beyond  the  scope  of  this 
Handbook. 


An  extensive  table  of  conical  flow  functions  is  given  in  the 
Princeton  Series,  Vol.  VII,  pp,  156-165  (Ref.  37). 
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702  Calculat ion  of  W ing  Characteristics 

The  lift,  moment,  and  drag  curves  of  representative  wing  plan- 
forms  are  presented  in  this  section  of  the  Handbook.  These  results  are 
based  on  the  method  of  supersonic  source  distributions  treated  in  Sub¬ 
sections  701  through  701.31  above.  When  the  wing  is  taken  at  an  angle 
of  attack,  the  flow  must  satisfy  the  Kutta  condition  (i.e.,  finite 
velocity  at  the  trailing  edge) .  A  trailing  edge  for  which  the  normal 
component  of  flow  is  supersonic  automatically  satisfies  this  condition. 
General  solutions  for  the  subsonic  trailing  edge  condition  have  not 
been  obtained.  Consequently,  the  lift  and  moment  data  presented  in 
this  section  are  restricted  to  supersonic  trailing  edges.  The  drag 
results  do  not  experience  this  limitation  because  they  are  restricted 
here  to  zero  angle  of  attack. 


702.1  Lift  and  Moment  Characteristics 

The  limitations  of  the  small  perturbation  linearized  theory 
for  supersonic  flow  restricts  the  applications  to  thin  wings  at  small 
angles  of  attack.  The  slopes  of  the  lift  and  moment  curves  with  angle 
of  attack  depend  on  the  wing  pi an form  only  and  not  on  the  camber  and 
thickness.  The  wing,  therefore,  can  be  treated  as  a  flat  plate,  and 
lift  (under  the  small  angle  of  attack  limitation)  is  assumed  equal  to 
the  normal  pressure  force  obtained  by  the  integration  of  the  pressure 
over  the  surface  of  the  wing.  Therefore,  the  lift  is  given  by 

L  =  |  (p^  -  pu)dS  , 

S 

where  the  subscripts  l  and  u  refer  to  the  lower  and  upper  surfaces 
respectively  and  S  is  the  wing  planform  area.  This  expression  for  the 
lift  is  rewritten  by  use  of  Eqs.  701.32-1  and  701.32-2  to  obtain  a  lift 
coefficient  expressed  in  terms  of  the  perturbation  velocity,  u: 


S 


The  ratio  of  the  perturbation  velocity  to  the  free  stream  velocity,  y, 

is  shown  by  Eqs.  701.42-1,  701.42-2,  and  701.42-3  to  be  proportional 
to  the  slope  of  the  wing  surface,  u>x.  The  upper  and  lower  surfaces  of 

a  flat  plate  wing  are  parallel  and,  under  the  small  angle  of  attack 
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limitation,  the  slopes  are  equal  to  (- a )  and  (a)  respectively,  where 
a  is  the  angle  of  attack  in  radians.  Therefore,  the  slope  of  the  lift 
curve  is  given  by 


(702.1-1) 


The  integral  expression  for  the  wing  pitching-moment-curve 

slope  is 


y 


i 

s 


(702.1-2) 


where  x  is  the  distance  from  the  center  of  the  element  of  area  to 

cp 

the  moment  axis,  and  cr  is  the  wing  chord  at  its  root  section. 


The  evaluations  of  Eqs.  702.1-1  and  702.1-2  are  accomplished 
most  conveniently  by  use  of  the  superposition  principle,  since  the  per¬ 
turbation  velocity,  u,  is  not  known  as  a  single  function  of  the  wing 
surface  coordinates.  The  triangular  wing  is  used  as  a  basic  configura¬ 
tion  for  which  basic  lift  and  moment-curve  slopes  are  obtained  by  the 
use  of  Eqs.  701.42-1  and  701.42-2,  or  701.42-3  in  Eqs.  702.1-1  and 
702.1-2  depending  upon  whether  the  leading  edge  is  supersonic  or  sub¬ 
sonic.  The  integrations  are  simplified  by  using  a  triangular  element 
of  area  with  its  apex  at  the  wing  apex  and  bounded  by  radial  lines  and 
a  segment  of  the  wing  boundary.  This  reduces  the  double  integrals  in 
Eqs.  702.1-1  and  702.1-2  to  single  integrals,  since  the  function  for 

remains  constant  along  a  radial  line.  The  new  variable  is  as 

Va  x 

used  by  Cohen  (Ref.  10). 


The  methods  described  in  the  preceding  paragraphs  have  been 
applied  by  various  researchers  to  the  determination  of  the  lift  and 
moment-curve  slopes  for  commonly  used  finite  wings  moving  at  supersonic 
speeds.  Formulas  resulting  from  these  analyses  are  presented  in  Cohen 
(Ref.  10)  and  Piland  (Ref.  11)  along  with  references  to  the  original 
sources.  Additional  applications  have  been  made  at  the  Defense  Research 
Laboratory  of  The  University  of  Texas  (Ref.  12).  This  latter  reference, 
which  also  includes  charts  presenting  the  lift  and  moment-curve  data 
from  all  of  the  above  sources,  is  used  as  the  main  source  of  the  numeri¬ 
cal  data  presented  in  this  section  of  the  Handbook. 
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702.11  Parameters  for  Lift  and  Moment  Data 

A  review  of  the  existing  data  and  of  the  methods  of  presenta¬ 
tion  has  led  to  the  choice  of  the  following  parameters  for  the  presen¬ 
tation  of  the  lift  and  moment  data: 


0  A 

0  cot  A 
*  % 


Mach  number-aspect  ratio  parameter 

Mach  number-leading  edge  sweepback  parameter 

Mach  number-lift-curve  slope  parameter 

Mach  number-moment-curve  slope  parameter 


The  variations  of  the  Mach  number-planf orm  parameters  0  A  and 
/3cot  A  with  Mach  number  for  certain  constant  values  of  A  and  A  are 
shown  in  Fig.  702.11-1.  The  planforms  of  swept  wings  with  unraked  tips 
can  usually  be  specified  by  three  parameters:  (1)  aspect  ratio,  A, 

(2)  leading  edge  sweepback  angle,  A  ,  and  (3)  taper  ratio,  A.  Commonly 
used  wing  planform  categories  are  defined  by  the  following  relations. 


A 


1  ~  A 
1  +  A 


cot  A 


(702.11-1) 


where  >  ,  =,  and  <  specify  rweptback,  unswept,  and  sweptforward  trail¬ 
ing  edges  respectively.  Planforms  having  fore  and  aft  symmetry  are 
defined  by  the  relation 

A  =  2  cot  A  (702.11-2) 

Illustrations  of  typical  configurations  so  defined  by  this  relation  are 
shown  in  Fig.  702.11-1  for  taper  ratios  of  A  =  0  and  A  =  0.5. 

It  has  been  shown  earlier  (see  figure  given  in  Subsection 
701.41)  that  regions  of  subsonic  and  supersonic  leading  edges  are  de¬ 
fined  by 

0  cot  A  =  1  ,  (702.11-3) 

where  <  ,  =  ,  and  >  specify  subsonic,  sonic,  and  supersonic  leading 
edge  respectively.  Regions  of  subsonic  and  supersonic  trailing  edges 
can  be  defined  by  simple  relations  between  the  planform  parameters  as 
given  on  the  following  page. 
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fi  cot  A 


Sweptback  trailing  edge: 
£  P  A 


>  4  ttt  +  e* 


4  Scot  A  < 

»r  l'-  /Scot  A"  >  **  • 


(702.11-4) 


where  the  signs  <  ,  =,  and  >  designate  subsonic,  sonic,  and  supersonic 
trailing  edge  conditions  respectively. 

Sweptforward  trailing  edge: 

flA  <  <  4  l"T"i  #  cot  A 

i  -  x  _  iScoli'orJi  ;  r~+  g  coTa - ’ 


t-rx -P* 


(702.11-5) 


where  <  ,  =,  and  >  specify  subsonic,  sonic,  and  supersonic  trailing 

edge  conditions  respectively.  Subsonic  and  supersonic  trailing  edge 
conditions  are  shown  in  the  following  figure. 


(a)  Subsonic,  Sweptback 
Trailing  Edge,  A  =  0 


(b)  Supersonic ,  Sweptforward 
Trailing  Edge,  A  =0.5 
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Wing  planforms  with  taper  ratios,  A  =  0,  0.5,  and  1.0  corres¬ 
ponding  to  various  combinations  of  $A  and  /Scot A  at  an  arbitrary  Mach 
number  of  M  =  2.24  (0=2)  are  shown  in  Figs.  702.11-2,  702.11-3,  and 
702.11-4  respectively  in  order  to  illustrate  the  variations  in  planforms 
associated  with  changes  in  the  planform  parameters.  The  short-dashed 
lines  shown  on  each  diagram  indicate  the  Mach  lines  for  the  Mach  number, 
M  =  2.24,  given  above.  The  apex  of  each  wing  leading  edge  designates 
the  0A  and  £cotA  coordinates  associated  with  that  wing  planform. 
Equations  702.11-1  and  702.11-2  show  that  the  unswept  trailing  edge  and 
symmetrical  planforms  are  associated  with  values  of  0  A  and  0  cot  A 
lying  along  a  radial  line  defined  by  constant  values  of  the  ratio, 


The  sonic  leading  edge  and  trailing  edge  boundaries  defined 


/SA 

£cot A  * 

by  Eqs.  702.11-3,  702.11-4,  and  702.11-5  are  represented  by  solid  lines 
shown  on  Figs.  702.11-2,  702.11-3,  and  702,11-4.  A  similar  figure  for 
trapezoidal  wings,  for  which  the  planform  parameters  are  0A  and  0tan  6 
is  shown  in  Fig.  702.11-5. 


702.2  Drag  Characteristics 

The  drag  of  finite  wings  is  comprised  of  three  parts:  namely, 
drag  due  to  thickness  of  the  airfoil  section  (usually  referred  to  as 
pressure  drag  or  wave  drag  at  zero  lift),  drag  due  to  lift  (or  induced 
drag) ,  and  skin  friction  drag. 

702.21  Discussion  of  Drag  Components 

The  pressure  drag  is  the  integral,  taken  over  the  wing  sur¬ 
faces,  of  the  drag  component  of  the  pressure  acting  on  those  surfaces. 
This  drag  is  known  to  be  zero  for  a  body  moving  at  subsonic  speeds  in 
a  i.  mviscous  incompressible  fluid.  In  an  actual  fluid,  however,  there 
exists  a  wake  due  to  the  viscous  boundary  layer  and,  due  to  the  pres¬ 
ence  of  this  wake,  the  full  pressure  recovery  over  the  after  portion 
of  the  wing  surface  is  not  realized.  Thus,  in  an  actual  fluid,  some 
pressure  drag  is  always  present.  Nevertheless,  efficient  wings  moving 
at  subsonic  speeds  produce  such  very  small  wakes  that  the  actual  pres¬ 
sure  drag  is  small  and  therefore  is  included  with  the  experimentally- 
determined  skin  friction  drag  to  make  up  the  well-known  profile  drag 
of  wings  at  subsonic  speeds. 

The  integrated  pressure  drag  for  finite  wings  moving  at  super¬ 
sonic  speeds  in  a  nonviscous  compressible  fluid,  however,  is  not  zero. 
This  component  of  the  drag  is  commonly  referred  to  as  the  pressure  drag 
(or  wave  drag)  at  supersonic  speeds  and  is  treated  separately  from  skin 
friction  drag.  Pressure  drag  at  zero  lift  is  presented  in  this  section 
of  the  Handbook. 


The  drag  due  to  lift  is  generally  referred  to  as  induced 
drag.  Induced  drag  is  made  up  of  two  parts  in  supersonic  flow.  One 
part  is  the  loss  of  momentum  in  the  general  flow  direction  caused  by 
the  deflection  of  the  airstream,  i.e.,  the  reaction  to  the  lift  force 
on  the  wing.  This  component  of  the  drag  in  supersonic  flow  corresponds 

(continued  on  following  page) 
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The  expressions  for  CD,  with  ^  given  by 

Eqs.  701.42-1,  701.42-2,  and  701 .42-3,  may  be  in¬ 
tegrated  immediately  to  yield  the  drag  of  a  wedge 
shaped  triangular  wing  with  an  unswept,  blunt 
trailing  edge  (see  Puckett,  Ref.  2).  However,  by 
the  principle  of  superporition  of  elementary  solu¬ 
tions  given  by  the  linearized  theory,  the  results 
given  above  may  be  extended  to  treat  the  double 
wedge  sections  discussed  in  Subsection  702.23.  For 

example,  the  drag  due  to  the  wing  segment  atmt'  of 
the  inset  figure  can  be  found  by  superimposing  the 

solutions  for  the  triangular  wing  att'  of  slope  u>x 

and  the  triangular  wing  mtt'  with  slope  -u)x.  The  Q 

drag  due  to  a  single  wing  segment  atmt’  (denoted 

hereafter  with  the  subscript  am)  is  given  as  a* 

function  of  A  and  A  of  the  segment  in  a  con- 
am  am 

venient  set  of  charts  in  Multhopp  and  Winter 
(Ref.  36)  and  reproduced  in  Fig.  702.22-1  of  this 
Handbook  by  courtesy  of  the  Controller,  Her 
Brittanic  Majesty's  Stationery  Office,  England.  By 
suitably  combining  solutions,  the  complete  double 
wedge  can  be  constructed.  If  the  functions  A Fam, 

AF  ,  and  AF  correspond  to  solutions  for  the 
mr  ar 

appropriate  segments,  the  total  drag  is  given  by 


Root  Chord 


eCD  .  <Xa  -  V 


(i)  (X  -  X  )  (X  -  X  ) 

c  'a  nr  '  m  r' 


(  AF  +  AF 
'  am  mr 


-  AFar>- 


These  results  are  easily  generalized  to  any  polygonal  cross  section. 


The  results  of  Puckett  for  the  triangular  planform  with  an 
unswept  trailing  edge  was  first  extended  by  Puckett  and  Stewart 
(Ref.  13)  to  wings  with  sweptback  and  sweptforward  trailing  edges. 
Similar  source  solutions  using  a  somewhat  different  approach  were  found 
by  Jones  (Ref.  14).  Margolis  (Refs.  15,  16,  and  17)  applied  the  method 
of  Jones  (Ref,  14)  to  the  evaluation  of  zero  lift  wave  drag  of  swept 
tapered  and  untapered  wings  having  a  double  wedge  profile.  Harmon  and 
Swanson  (Ref.  18)  and  Harmon  (Ref.  19)  applied  the  same  methods  to  the 
evaluation  of  the  wave  drag  for  wings  having  a  symmetrical  biconvex 
profile.  Chang  (Ref.  20)  applied  the  Von  Karman  integral  method  to  the 
analysis  of  the  flow  about  a  general  swept  tapered  planform  wing  at 
supersonic  speed.  The  wave  drag  expressions,  derived  by  these  writers 
and  used  for  the  computation  of  drag  values,  are  long  and  complicated. 
These  formulae  are  published  in  available  sources  and  are  not,  there¬ 
fore,  reproduced  in  this  Handbook. 


Reference  is  made  here  to  a  comprehensive  survey  by  Lawrence 
(Ref.  21)  of  sources  of  supersonic  drag  theory.  Drag  data  obtained 
from  the  original  sources,  supplemented  by  the  results  of  additional 
computations,  are  presented  in  Lawrence  (Ref.  21)  in  chart  form  for  a 
wide  range  in  planforms  of  double  wedge  and  biconvex  profile  wings. 
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702.23  Discussion  of  the  Curves  for  Lif t .  Moment ,  and  Zero  Lift  Drag 


Values  of  the  lift-curve  and  the  moment-curve  slopes  for  a 
range  of  several  commonly  used  finite  wing  planforms  in  supersonic  flow 
have  been  computed  at  the  Defense  Research  Laboratory  of  The  University 
of  Texas.  These  data  are  presented  in  the  form  of  a  modified  lift- 
curve  slope,  yS  C.  ,  and  a  modified  moment-curve  slope,  /3  C  ,  plotted 

a 

against  y3  cot  A  for  certain  designated  constant  values  of  /3A.  The 
charts  presenting  values  of  the  modified  lift-curve  slope,  /9  C,  ,  are 

a 

extended  into  the  subsonic  sweptforward  trailing  edge  region  by  means 

of  the  reversibility  theorem  (see  discussion  in  Subsection  702.3). 

Since  this  theorem  does  not  apply  to  the  moment  characteristics,  there 

is  no  data  available  for  the  moment-curve  slope,  BC  ,  in  the  subsonic 

m 

a 

sweptforward  trailing  edge  region.  Lift-curve  and  moment-curve  slopes 
for  the  subsonic  sweptback  trailing  edge  condition  have  not  been  ob¬ 
tained  (see  Subsection  702). 


The  zero  lift  pressure  drag  data  for  these  wing  planforms 
with  symmetrical  double  wedge  profiles  are  presented  in  this  Handbook 
in  the  form  of  a  Mach  number-drag-thickness  parameter,  y3C„  /  .  2.  This 


parameter,  / .  2,  is  plotted  versus  (3  cot  A  for  constant  values  of 

V(b 


j3 A  as  was  done  for  the  CT  and  C  data. 
H  La  ma 


The  lift-curve  slopes,  the  moment-curve  slopes,  and  the  data 
for  pressure  drag  at  zero  lift  for  finite  planforms  where  A  =  0,  0.5, 
1.0,  and  trapezoidal  wings  are  presented  in  Figs.  702.23-1  through 

702.23- 16.  Similar  data  for  rectangular  wings  can  be  found  in  Figs. 

702.23- 15  and  702.23-16  (tan  6  =  0),  and  Fig.  702.23-12  (  £  cot  A  =  °o). 
Composite  diagrams  showing  the  areas  of  solution  for  the  lift-curve 
and  moment-curve  slope  curves  and  drag  curves  of  these  planforms  are 
included  in  this  group  of  figures  (Figs.  702.23-1,  702.23-5,  702.23-9, 
and  702.23-14).  Primary  reference  sources  for  the  data  are  cited  in 
these  figures  for  the  convenience  of  the  reader. 


702.3  Reversibility  Theorem 

The  assumptions  of  the  linearized,  potential-flow  theory  make 
possible  the  statement  of  a  simple  reversible  flow  theorem.  Several 
authors  have  treated  this  problem  under  various  restrictions  (see  Refs. 
22  through  26).  Brown,  (Ref.  27)  has  shown  that  the  lift-curve  slope 
and  the  thickness  drag  remain  the  same  when  any  airfoil,  or  system  of 
airfoils,  is  reversed  so  long  as  the  Kutta  condition  of  finite  velocity 
at  the  wing  trailing  edge  applies.  This  theorem  has  enabled  the  exten¬ 
sion  of  the  data  presented  in  this  Handbook  into  certain  areas  not 
otherwise  calculable.  It  allows  one  to  determine  the  characteristics 
for  those  planforms  which  can  be  identified  with  an  equivalent  configu¬ 
ration  for  which  the  lift-curve  slope  and  the  thickness  data  are  avail¬ 
able  in  the  Handbook.  Examples  illustrating  these  applications  are 
included  in  the  following  subsections. 
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703  Numerical  Examples 

The  evaluation  of  finite  wing  characteristics  by  the  use  of 
the  data  presented  in  this  section  of  the  Handbook  is  illustrated  in 
the  following  examples. 


703.1 


Example  1 :  Wing  with  a  Supersonic  Leading  and  Trailing  Edge 


The  determination  of  the  theoretical  lift-curve  slope,  the 
moment-curve  slope,  and  the  pressure  drag  coefficient  for  the  following 
finite  symmetrical  double  wedge  wing  is  illustrated  in  this  example. 

Given : 


Mach  number  M  =  3.16 

Wing  specifications 

Leading  edge  sweep  angle  A  =  60° 


Thickness  ratio 


Tip  chord 


Wing  span 


t  =  °‘09 


ct  '  2  cr 


b  =  f  cr 


Computations : 


-it;t 


0  = yvf  -  i 

-  3.000 

*  cot  A  =  0.577 

4  b 

A  =  g-  —  (given  from  the  geometry  of 
r  the  inset  figure) 

=  2 

Thus,  /3A  =  6.000 

and  ficot  A  =  1.731 

Figure  702.11-3  indicates  supersonic  leading  and  trailing  edges. 

0C.  =  4.11  (from  Fig.  702.23-6) 

a 

(3  C  =  -4.03  (from  Fig.  702.23-7) 


=  4.59  (from  Fig.  702.23-8) 


Results ; 


C,  =  1.37  per  radian  (  0.0239  per  degree) 

Cm  =  -1.34  per  radian  (-0.0234  per  degree) 
a 

CD  = 


0.0124 
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703.2  Example  2:  Wing  with  a  Supersonic  Leading  Edge  and  a  Subsonic 
Trailing  Edge 


The  determination  of  the  theoretical  lift-curve  slope,  the 
moment-curve  slope,  and  the  wave  drag  coefficient  for  a  finite  symmetri¬ 
cal  double  wedge  wing  is  illustrated  in  this  example. 


Given : 


Mach  number  M  =  1.25 

Wing  specifications 

Leading  edge  sweep  angle  A  =  30° 

Thickness  ratio  ~  =  0.09 

Triangular  panel  A  =  0 


Wing  span 


Computat:  ns: 


0  = 


(3  cot  A  = 
A  = 


Thus,  (3  A  = 
and  (3  cot  A  = 


|/m2  -  1 

0.750 

1.732 

ty  i. 

—  (given  from  the  geometry  of 
r  the  inset  figure) 

2 

1.500 

1.300 


r 


Reference  to  Fig.  702.11-2  shows  that  this  wing  lies  in  the  subsonic 
trailing  edge  region.  Figure  702.23-1  indicates  that  the  lift-curve 
slope  and  wave-drag  coefficient  can  be  computed  by  the  application  of 
the  reversibility  theorem  referred  to  in  Subsection  702.3. 

This  theorem  simply  states  that  the  lift-curve  slope  and  wave- 
drag  coefficient  given  by  the  linearized  theory  for  a  given  wing  remain 
the  same  when  the  flow  direction  is  rotated  through  180  degrees,  pro¬ 
vided  the  Kutta  condition  at  the  trailing  edge  is  satisfied.  This  con¬ 
dition  is  automatically  satisfied  for  wings  at  zero  lift  as  is  true  for 
the  wing  data  presented  herein. 

The  specified  leading  edge  sweep  angle  of  the  wing  in 
Example  2  is  A  =  30°  and  A  =  0.  The  trailing  edge  sweep  angle  can  be 
seen  to  be  given  as 


cot  A  T  =  -0.703  . 

The  negative  sign  indicates  a  sweptforward  trailing  edge  for  the  speci¬ 
fied  wing  and  therefore  the  reverse  flow  wing  would  have  a  sweptback 
leading  edge.  The  computations  for  the  reverse  flow  wing  are  made  in 
the  following  manner. 
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703.21  Use  of  the  Reversibility  Theorem  to  Extend  the  Curves 


Reverse  flow  wing  computations: 
jS  cot  A  =  0.527 
P  A  =  1.500 

=  2.41  (from  Fig.  702.23-2) 


0C. 


a 


I3C 


m 


a 


^CD 

t  2 

<£> 

CL 

a 

CD 


-  reversibility  theorem  not  applicable 
=  3.74  (from  Fig.  702.23-4) 

=  3.21  per  radian  (0.054  per  degree) 

-  0.0404 


The  curves  in  Figs.  702.23-2  and  702.23-4  have  been  extended 
into  the  subsonic  trailing  edge  region  by  the  method  illustrated  in 
this  example. 


703.3  Example  3^:  Wing  with  £  Sweptf orward  Leading  Edge 


The  determination  of  the  theoretical  lift-curve  slope  and  the 
wave-drag  coefficient  for  wings  with  a  sweptforward  leading  edge  is 
illustrated  in  this  example.  These  data  are  obtained  by  computations 
for  the  equivalent  reverse  flow  wing  as  demonstrated  by  the  reversibility 
theorem : 


Given : 

Mach  number  M 

Wing  specifications 

Leading  edge  sweep  angle  A 

Thickness  ratio  t^ 

c 

Tip  chord  ct 

Wing  span  b 


1.80 
-45° 
0.09  V 


\  cr 

hr 


r 


\ 


/■ 

\ 


Computations : 

cot  A  =  -1.00 

From  the  given  values  of  A  and  X, 
cot  A  =  -0.60 


489030  0  -59  -3 
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Since  cot At  is  negative  (thereby  indicating  a  sweptforward 

trailing  edge)  the  equivalent  reverse  flow  wing  would  have  a  sweptback 
leading  edge.  The  computations  for  the  equivalent  reverse  flow  wing 
are  made  in  the  following  subsection. 


703.31  Reverse  Flow  Wing 


Computat ions : 

(3  =  a/m2  -  1 


= 

1.500 

A  = 

4  b 

3  cr 

(given  from 
of  the 

the  geometry 
inset  figure) 

= 

2 

Thus,  /3  A  = 

3.00 

(3  cot  A  = 

0.900 

and  /3C,  = 

a 

3.47 

( from  Fig. 

702.23-6) 

0CD 
♦  2 

5.35 

(from  Fig. 

702.23-8) 

✓ 


Results : 

C.  =  2.31  per  radian  (0.0403  per  degree) 

La 


C 


D 


0.0289 
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704  Comparison  of  Theory  and  Experiment 

A  systematic  comparison  of  the  measured  and  the  theoret ic  .  1 i 
calculated  aerodynamic  characteristics  of  a  large  number  of  finite  :ags 
at  supersonic  speeds  has  been  made  by  Vincenti  at  the  Ames  Aeronautical 
Laboratory  of  the  National  Advisory  Committee  for  Aeronautics  ( Rel .  28). 
Approximately  thirty  wing  models  were  tested  at  a  Mach  number  of  1.53 
in  the  experimental  investigation  in  which  the  lift,  pitching-moment, 
and  drag  characteristics  were  determined.  Figures  704-1  through  704-6 
in  this  section  of  the  Handbook  present  a  graphical  summary  of  these 
results  based  on  Figs.  4  through  9  of  the  reference  cited  above.  The 
comparisons  given  in  these  figures  should  give  a  convenient  basis  for 
estimating  the  accuracy  and  reliability  of  engineering  predictions 
from  the  linearized  theory. 

Comparisons  of  lift-curve  slopes  as  determined  by  measurement 
and  theory  for  variations  in  aspect  ratio  and  sweep  angle  are  shown  in 
Figs.  704-1  and  704-2.  The  wings  in  these  instances  had  isosceles  tri¬ 
angle  sections  of  5  per  cent  maximum  thickness.  The  agreement  with 
respect  to  aspect  ratio  is  so  good  that  Vincenti  suggests  that  the  ef¬ 
fects  of  viscosity  and  support-interference  must  be  almost  completely 
compensating  for  the  wings  tested.  The  trends  in  the  effects  of  sweep 
on  the  lift-curve  slopes  (Fig.  704-2)  as  predicted  by  the  linear  theory 
are,  therefore,  confirmed  by  the  experimental  results.  The  symmetry  of 
the  lift-curve  slope  with  variation  in  sweep  angle  verifies  the  theoreti¬ 
cally  predicted  invariance  of  the  lift-curve  slope  with  the  flow  direc¬ 
tion  reversed.  However,  the  generally  acceptable  agreement  in  over-all 
lift  characteristics  must  not  be  construed  to  imply  satisfactory  agree¬ 
ment  with  regard  to  load  distribution. 

Similar  comparisons  for  the  moment-curve  slope  are  shown  in 
Figs.  704-3  and  704-4.  The  agreement  in  these  instances  is  neither 
qualitatively  nor  quantitatively  good.  The  lack  of  symmetry  with  re¬ 
spect  to  sweepback  and  sweepforward  (Fig.  704-4)  verifies  the  non¬ 
reversibility  prediction  for  moment  characteristics.  Vincenti  suggests 
that  the  development  of  second-order  wing  theory  may  be  necessary  for 
solution  of  the  pitching-moment  prediction  problem. 

The  effect  of  variation  in  sweep  on  the  minimum  drag  is  shown 
in  Fig.  704-5.  Again  the  symmetry  verifies  the  conclusions  of  reversi¬ 
bility  theorems.  The  amount  of  the  increase  in  drag  with  increase  in 
sweep  predicted  by  linear  theory  is  not  observed  experimentally.  The 
difference  between  the  measured  and  theoretical  drag  at  low  sweep  angles 
is  consistent  with  a  reasonable  estimation  of  the  friction  drag  which 
must  be  added  to  the  theoretical  wave  drag  to  obtain  the  total  drag. 
However,  it  appears  on  the  same  basis  that  linear  theory  predicts  too 
high  a  drag  as  the  wing  sweep  approaches  coincidence  with  the  Mach  line 
direction.  Theoretical  results  for  the  subsonic  trailing  edge  condition 
are  not  available,  but  the  generally  predicted  reduction  in  drag  as  the 
wing  is  swept  behind  the  Mach  cone  is  verified  by  the  experimental  re¬ 
sults.  Figure  704-6  shows  a  comparison  for  results  on  double  wedge 
profile  triangular  wings  of  varying  maximum  thickness  position.  The 
addition  of  an  estimated  skin  friction  to  the  wave  drag  produced  a  rea¬ 
sonably  satisfactory  agreement. 

It  must  be  concluded,  therefore,  that  linearized  theory  data 
can  be  used  only  as  a  general  guide  for  design  predictions. 
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705  Limitat ions  of  the  Linearized  Theory 

The  effect  of  viscosity,  the  effect  of  wing-body  interference, 
and  higher  order  corrections  to  the  linearized  theory  must  be  accounted 
for  in  any  design  applications.  The  reader  is  referred  to  Hilton 
(Ref.  29)  and  Shapiro  (Ref.  30)  for  a  more  general  discussion  of  the 
three-dimensional  wing  problem  at  supersonic  speeds. 


705.1  The  Effect  of  Viscosity 

The  drag  data  presented  in  this  section  do  not  include  any 
viscous  effects.  In  the  absence  of  reliable  experimental  information 
on  wing  drag  at  supersonic  speeds,  a  useful  estimation  can  be  arrived 
at  by  the  addition  of  a  calculated  average  skin  friction  coefficient 
such  as  presented  in  Section  6  of  the  Handbook  (Two-Dimensional  Air¬ 
foils)  . 


705.2  Wing-Body  Interference 

An  important  factor  which  is  always  present  to  modify  the  pre¬ 
dicted  summation  of  characteristics  from  the  component  parts  is  wing- 
body  interference  effects.  In  general,  the  evaluation  of  these  effects 
must  be  obtained  through  precise  experimental  and  analytical  studies. 

The  importance  and  extensiveness  of  this  subject  warrants  its  presenta¬ 
tion  in  a  separate  section  of  the  Handbook  (Section  9:  Mutual  Inter¬ 
ference  Phenomena)  where  this  problem  is  discussed  in  detail. 


705.3  Higher  Order  Theories 

The  need  for  second  order  theory  for  finite  wings  at  super¬ 
sonic  speeds  was  suggested  earlier  as  a  result  of  the  unsatisfactory 
prediction  of  pitching-moment  (see  Subsection  704  and  Figs.  704-3  and 
704-4) .  The  usefulness  of  extending  the  first  order  theory  to  a  higher 
order  theory  is  somewhat  questionable  because  of  the  possible  presence 
of  large  viscous  effects.  The  magnitude  of  the  effects  of  higher  order 
approximations  in  two-dimensional  flows  is  thoroughly  treated  in  Section 
6  of  the  Handbook  (Two-Dimensional  Airfoils). 

Van  Dyke  (Ref.  31)  discusses  the  development  of  a  second  order 
theory  in  some  detail,  and  Lighthill  (Ref.  35)  presents  a  general  dis¬ 
cussion  of  higher  order  approximations  for  supersonic  three-dimensional 
wing  theory.  However,  results  of  immediate  interest,  beyond  those  of 
the  linearized  theory,  do  not  appear  to  be  available. 
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702*22  Figure  1  Chart  for  the  function  A  F  to  be  used  in  the  calculation  of 
pressure  drag  of  X  =  0  planforras  (Taken  from  Ref.  36). 
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702.23  Figure  2  Generalized  curves  of  the  wing  lift-curve-slope  parameter 
for  tapered  wings  with  swept  leading  edges,  A  =  0. 
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702.11  Figure  3  Variations  in  planforra  associated  with  changes  in  the 
planform  parameters,  @A  and  /3cotA  ;  X  =  0.5. 
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702*23  Figure  6  Generalized  curves  of  the  wing  lift-curve-slope  parameter 
for  tapered  wings  with  swept  leading  edges,  X  -  0.5. 
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2 .23  Figure  7  Generalized  curves  of  the  wing  moment -curve-slope  parameter 
for  tapered  wings  with  swept  leading  edges,  X-  0*5. 
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702.23  Figure  8  Generalized  curves  of  the  wing  pressure-drag 
parameter  at  zero  lift  for  tapered  wings  with 
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swept  leading  edges;  double  wedge  profile,  maximum  thickness 
at  0.5  chord,  \  =  0.5. 


Sup«r»onic  L.E.and  T.E 


August  1957 


A  =  1.0  Planforms 


702.11  Fig.  4 


$  Cot  A 

702.11  Fig.  4  Variations  in  planform  associated  with  changes  in  the 
planform  parameters /9  A  and/3cotA  ;  X=  1-0 
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702.23  Fig.  9  Three-Dimensional  Airfoils  August  1957 
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702.23  Figure  10  Generalized  curves  of  the  wing  lift-curve-slope  parameter 
for  untapered  wings  with  swept  leading  edges,  X  -  1.0. 
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702.23  Figure  11  Generalized  curves  of  the  wing  moment-curve-slope  parameter 
for  untapered  wings  with  swept  leading  edges,  X  -  1.0. 
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702.23  Figure  13  Generalized  curves  of  the  wing  pressure-drag 
parameter  at  zero  lift  for  untapered  wings  with 
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=  1.0  Planforms  (Biconvex  Profile)  702.23  Fig.  13 
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swept  leading  edges;  biconvex  parabolic-arc  profile,  maximum 
thickness  at  0.5  chord,  X  =  1.0. 
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Trapezoidal  Planforms  702.23  Fig.  14 
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Figure  16  Generalized  curves  of  the  wing  moment-curve-slope 
parameter  for  trapezoidal  wings. 


NAVORD  REPORT  1488  (Volume  3) 

704  Fig,  1  Three-Dimensional  Airfoils  August  1957 


August  1957 


NAVORD  REPORT  1488  (Volume  3) 

704  Fig.  3  Three-Dimensional  Airfoils 
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Comparison  of  Theory  and  Experiment 
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Comparison  of  Theory  and  Experiment 


704  Fig.  6 


Position  of  Maximum  Thickness, 
Percent  of  Chord 


704  Figure  6  Effect  of  position  of  maximum  thickness  on  the  minimum 
drag  of  an  uncambered  triangular  wing;  double  wedge  profile, 

M  =  1.53,  -  =  0.05,  A.  =  0,  a  =  0,  and  R  =  750,000 
oo  c  e 

(based  on  the  mean  geometric  chord  of  the  wing).  (Ref.  28.) 
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